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Lattice gauge calculations (1) for static (heavy) quarks sup- 
port the notion that the interquark potential in quantum chro- 
modynamics (QCD) behaves as V (r) - \r for large r. Indeed, 
the linear potential has long been used in phenomenological 
nonrelativistic quark models of baryons and mesons (2, 3). 
Meson spectroscopy in particular is successfully described by 
a linear potential at large r, modified by spin- and colour- 
dependent Coulomb forces at small r. Most calculations with 
the linear potential are earned out in coordinate space. This is 
the simplest procedure for heavy -quark systems, which can per- 
haps be considered as nonrelativistic; however for light-quark 
systems it would be desirable to have a relativistic treatment. 
Bound-state equations in relativistic systems (4) are generally 
much easier to solve in momentum space, and thus we are led 
to consider, as a starting point for the relativistic case, the 
Schrodinger equation for two scalar particles interacting by a 
linear potential. The methods developed will generalize rela- 
tively straightforwardly to relativistic treatments. 

To summarize: here, we treat the Schrodinger equation for 
a linear r-space potential. The method is for the most part 
straightforward, the only difficulty arising from the singularity 
of the kernel at the origin of momentum space. Previous treat- 
ments (5) have usually been approximate in the sense that the 
singularity was handled by screening the r-space pontendal: 

[1] V{r)~\rt~ v 

What has perhaps not been generally appreciated is that the 
limit t| — » 0 can be taken analytically. Previous treatments keep 
the parameter t\ finite, leading to some uncertainty as to the 
nature of the calculated eigen values and wave functions. In this 
connection, recall that the screened linear potential does not, 
strictly speaking, possess true bound states, instead it has scat- 
tering resonances, which for low energy approximate the bound 


states of the unscreened potential. We will extract the limitof 
zero screening analytically, using a subtraction technique. The 
resulting subtracted integral equation is relatively easy to han- 
dle numerically. An alternative procedure, not employing any 
subtraction, and leading to a different integrodifferential equa- 
tion is presented in ref. 6 Our approach is easy to implement 
and generalizes without difficulty to higher partial waves. The 
Schrodinger equation for the hh partial wave is (with the 
inhomogeneous term already omitted, as it will not contribute 
to the bound states in the limit of zero screening) 


[2] ^ 4»X p) + j Vfp.p’)4>{p')p' 2 dp' * EMP) 


Here p. * + m,) is the reduced mass and V,, giver 
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is the Ah partial-wave component of the Fourier transform o 
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The variable y is given by: 

v + £ t £ 


[5] 


2 PP' 


NORBURY ET At 




Q',(y) and Q",(y) are the first and second derivatives (with 
respect to y ) of the Legendre function of the second kind. To 
illustrate the method we specialize to r-waves, where we find 
by contour integration 

l6J / np.p , )p' 1 dp' « “ / dp' |&(y) + ^CS(y) J 



Note that when -q * 0, Qtfy) and Qtfy) have double and 
quadruple poles, respectively, at p' * p, so that their integrals 
do not exist separately. Nevertheless, the two terms added 
together produce a function with an integrate singularity. This 
is illustrated in Fig. 1 , which shows the kernel as a function of 
p' for fixed p. One observes that there is a central maximum 
atp' = p with height scaling as 1 /y| 2 , flanked by two minima 
atp' -p ± 2t) whose heights also scale with Urf. The integral 
vanishes [6] and this allows us to rewrite the Schrodinger equa- 
tion in subtracted form 


[7] 


L 

2H 


<t> 0 (P) + ~ 
IT 



+ V 


ggCy) 1 

(pp') 1 J 


X [<J>0 (/>') - ^p))p' 2 dp' = £<h 0 (p) 

The limit tj — ► 0 now exists, and may be extracted by splitting 
the region of integration to isolate the singularity. We write 

hi / ¥ ~ eiwl - *,(?)] 


r ♦ p - r 
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= i 4 + fl + C 

The limit sp ± 4tj are chosen so that ail three extrema of the 
kernel lie in the middle region B . The explicit forms of the 
Legendre functions are 

«o> - rT7 



Fig. 1 . The singularity structure of the kernel is shown for finite 
n — 0.075 with fixed p * 2. 
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and 
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It is clear that for p' * p, as is the case in the integrals A and 
C, the limit tj— ► 0 is innocuous, and may be taken immediately, 
indeed one has 

[9] lim [A + C] 

*1—0 
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where P denotes as usual, the Cauchy principal value of the 
integral, which has been made well defined by the subtraction. 
The term B must be handled with care, however, since p' * p 
inside the region of integration. Assuming 4>(p') is analytic in 
the neighborhood of p, and making an obvious change of 
variable we find 
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Scaling out 4q then results in 
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Table I. Energy eigenvalues in GeV for l » 0. m, » m, » 1.5 GeV, and X » 5 GeV 2 





N 



Exact 

8 

10 

12 

14 

16 

18 

E, 

5.973 

5.972 

5.972 

5.972 

5.972 

5.972 

5.971 

E t 

10 46* 

10.444 

10.443 

10.443 

10443 

10.443 

10.441 

E, 

14.389 

14.114 

14.111 

14.104 

14.104 

14.104 

14.101 

E t 

18.646 

17.452 

17.378 

17.341 

17.335 

17.335 

17.335 

E , 

23.402 

21.125 

20.397 

20.351 

20.294 

20.293 

20.291 

E t 

27.206 

25.683 

23.440 

23.281 

23.072 

23.053 

23.046 

E, 

33.032 

31.269 

27.274 

26.059 

25.842 

25.648 

25 646 

E, 

44.374 

36.224 

32.113 

29.032 

28.789 

27.947 

28.119 

E, 


40 519 

38.146 

33.051 

31.561 

30.194 

30.488 

E, o 


51.774 

45.309 

38.067 

34.428 

33.340 

32.769 

E„ 



49.940 

44.286 

38.517 

36.489 

34.972 

E tl 



58.588 

51.893 

43.615 

37.309 

37.109 


The contribution of the second term in B1 clearly vanishes since it is not singular atp' = p, the analysis of B 2 is similar, and we 
conclude that B tends to zero. Therefore the limiting form of the equation is 

1121 - 7?' f *{<^7?] * E * M 

We now discuss the numerical solution of [12], which is not yet a completely trivial matter, since care must be taken to obtain 
the Cauchy principal value. In this respect there is a difference between the linear potential and the Coulomb potential, the latter 
giving rise to a logarithmic singularity. For the Coulomb potential, the method used in the literature (7) is to write the Coulomb 
analog of [12] directly, for example, using Gaussian quadrature, as a matrix equation. Since the singularity is only logarithmic 
this method is successful for the Coulomb potential. Here, such an approach is not feasible. Instead, we expand 4> 0 in a suitable 
set of basis functions 

s 

[13] 4>„ (P) = 2 C^,(p) 


Inserting this expansion in [12], multiplying by p 2 g m (p) and integrating overp, we obtain 

[14] 2 C,jj ^gm(p)g,(p)4p + “ / gm(p)[g'(p') - &.(/>)] dp' dp = E 2 C* j P 2 gm(p)g,(p)ty 


which is just the matrix equation 

[i5] 2 - £ 2 cj: m 


The double integral over p and p' is performed by changing to 
variables ( p ' + p) and (p' - p). The singularity is in the 
integral over (p' — p), so this is carried out first using Gaussian 
quadrature with an even number of points. This type of inte- 
gration yields the Cauchy principal value automatically (8). A 
convenient set of functions g( p) is 


116] g m (p) 


I 

{n 2 /N) 2 + p 4 


where N is the maximum number of functions used in expansion 
[13]. Figure 2isa3D plot of the kernel of [14], showing clearly 
the cancellation that leads to the principal value. Using the 
above method, we have calculated both eigenvalues and eigen* 
vectors. In Table 1 the first 12 eigenvalues are listed. We used 
m, * * 1 .5 GeV and the string tension X * 5 GeV 2 . One 

can see that the lower eigenvalues converge nicely as the num- 
ber of functions is increased. We compare these with the eigen- 
values obtained from a coordinate space calculation (integrating 
the equation out from r«0 and in from large r, and matching 


t integrand 



Fig. 2. A three-dimensional figure of the subtracted, regulated inte- 
grand; ti - 0.075. The cancellation that produces the Cauchy prin- 
ciple value is evident. 

the logarithmic derivatives at the classical turning point), in 
Table 1. The calculated eigenfunctions also agree with the 
coordinate-space calculation. 
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In conclusion, we have treated the problem of two nonrela- 
tivistic, scalar panicles interacting via a linear potential in 
momentum space. The relevant Schrodinger equation has a sin- 
gular kernel. We have shown how after regulating the singu- 
larity by exponentially screening the r-space potential, the 
severity of the singularity can be reduced by a suitable sub- 
traction, and the limit of zero screening extracted analytically. 
To the best of our knowledge, this point has not been generally 
understood in the literature. The limiting form of the equation 
has been treated numerically, and the results are in good agree- 
ment with more straightforward coordinate space calculations. 
Relativistic equations involving linear potentials involve sim- 
ilar singularities, so that the methods developed here will be 
applicable. We intend to study the relativistic quark-antiquark 
problem in the future. The method presented here can be 
generalized to higher partial waves without undue difficulty. 
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An optical potential abrasion-ablation collision model is used to calculate hadronic dissociation 
cross sections for one, two, and three nucleon removal for the first time for a 14.6 A GeV 28 Si 
beam fragmenting in aluminum, tin, and lead targets. These estimates are compared with recent 
semi-inclusive measurements. Significant differences between some calculated and measured 
semi-inclusive cross sections exist which cannot be resolved without measurements of the exclusive 
channel hadronic cross sections. Calculations for each exclusive reaction channel contributing to 
the semi-inclusive cross sections are presented and discussed. 


-Recently, the E814 Collaboration at Brookhaven Na- 
tional Laboratory (BNL) made a very detailed experi- 
mental study of the breakup of silicon beams at relativistic 
energies (.Eiab 1 " 14.6/4 GeV or Tub “l 3.7/4 GeV) using 
the Alternating-Gradient Synchrotron. 1 They reported 
cross-section measurements of one, two, and three nucleon 
removal by aluminum, tin, and lead targets from both 
electromagnetic and hadronic dissociation processes. For 
the electromagnetic dissociation (EMD) process, mea- 
surements of individual exclusive channel contributions 
were reported. Comparisons of measurements for 1 p and 
In removal with calculated values obtained using the 
Wcizsacker-Williams method of virtual quanta 2 were 
made, and good agreement was obtained. More recently, 
Llope and Braun-Munzinger 3 extended the EMD analysis 
to include multiple excitations of the giant dipole reso- 
nance coupled with fragmentation probabilities obtained 
from the standard statistical model of nuclear decay. 
They then use this extended calculational framework to 
predict exclusive EMD cross sections for many of the 
channels measured by the E814 Collaboration. 

For the measured hadronic dissociation channels, how- 
ever, no detailed analyses have been reported. In Ref. I, 
simple comparisons between semi-inclusive measurements 
and a recent parametrization 4 of 1 p and \n geometrical 
calculation of single nucleon removal 5 were made. In this 
work, we analyze the hadronic dissociation of silicon pro- 
jectile nuclei by aluminum, tin, and lead targets using an 
optical potential abrasion-ablation collision model which 
includes contributions from frictional-spectator interac- 
tions. 6 This model is used to calculate exclusive cross sec- 


tions. Although no exclusive experimental hadronic cross 
sections were reported in Ref. 1 (the only exclusive cross 
sections reported were due to EMD), these calculated re- 
sults are presented to stimulate interest in their experi- 
mental measurement and to facilitate further discussion in 
the semi-inclusive cross-section analysis. 

The abrasion portion of this formalism was recently 
used to successfully describe single nucleon emission in 
relativistic nucleus-nucleus collisions. 7 Predictions of ha- 
dronic cross sections for the exclusive reaction channels 
measured in Ref. I are presented. Semi-inclusive cross 
sections, obtained by summing the appropriate exclusive 
channels, are presented and compared with the measured 
values reported in Ref. 1. Reasonably good agreement is 
obtained for the xp Ox “l, 2, 3) channels. However, for 
the yn (y “1, 2) channels, the agreement is not as good, 
with the calculations generally overestimating the experi- 
mental data. Comments concerning the difficulties in 
resolving these differences are made, and the need for ex- 
clusive measurements of these hadronic cross sections is 
pointed out. 

In the optical potential formalism, 6 the abrasion cross 
section for removal of m nucleons is 


CTabr(^PF) 


Ap 

m 


JVmi -P(b)] m lP(b)] Arf , (1) 


where 

P(b) " exp[ — /ircr(e)/(A)] , (2) 


with 
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In Eqs. (2) and (3), b is the impact parameter, e is the 
two-nucleon kinetic energy in their center-of-ma ss frame, 
z o is the target center-of-mass position in the projectile 
rest frame, denotes the target-nucleus internal coordi- 
nates, and y is the projectile-nucleon-target-nucleon rela- 
tive coordinate. Methods for obtaining the appropriate 
nuclear distributions p, (i—P 9 T) and constituent- 
averaged nucleon-nucleon cross sections a(e) are given in 
Ref. 8. Values for the diffractive nucleon-nucleon scatter- 
ing slope parameter B(e) are obtained from the parame- 
trization in Ref. 9. The Pauli correlation correction de- 
rived in Ref. 8 is neglected here because it is negligible for 
the peripheral collisions 6 being considered in this work. 

Since the abraded nucleons consist of protons and neu- 
trons, a prescription for calculating the prefragment 
charge dispersion is needed. The three available choices 
are completely correlated, 10 hypergeometric (completely 
uncorrclated), 11 and a model based upon the zero-point 
vibrations of the giant dipole resonance. 12 For the present 
work, we have chosen to implement the hypergeometric 
model, which assumes that there is no correlation at all 
between the neutron and proton distributions. For few 
nucleon removal processes, such as are being investigated 
here, the calculated results are not particularly sensitive to 
any of these particular charge dispersion methods. 13 For 
example, all three methods yield identical charge disper- 
sion results for single nucleon abrasions from self- 
conjugate nuclei. If z out of the original Z projectile nu- 
cleus protons are abraded along with n out of the original 
/V projectile neutrons, then the abrasion cross section be- 
comes 


<*zbr(ZpF,Apf) 


where 



A p 


m 


(4) 


where P c K is the BCV probability that an abraded nu- 
cleon escapes without undergoing any frictional-spectator 
interactions. 5 For the reactions considered herein, 
P«c~ 0.7. 

Depending upon the excitation energy, the excited 
prefragment will decay by emitting one or more nucleons, 
composites, or gamma rays. The probability a, ; (p) for 
formation of a specific final fragment of type i as a result 
of the deexcitation of a prefragment of type j which has 
undergone p frictional-spectator interactions is obtained 
using the EVA-3 computer code. 12 For m*l or 2 and 
p “0 (no FSI), the values of E cxc are less than 3 MeV for 
all targets and no particle emission occurs. Hence, the 
calculated cross sections for 27 Si + n, 26 Si + 2n, 27 Al+p, 
and 26 Mg+2/7 arise solely from the abrasion process. 
Whenever one or more FSFs (p“l, 2) occur for a frag- 
menting silicon nucleus, an additional (average) excita- 
tion energy of 31 MeV per FSI (computed using the 
methods of Ref. 11) is deposited in the prefragment. 
When these resultant excitation energies are used as in- 
puts into the EVA-3 code, the cross sections for the 
25 Mg + 2 p+n and the 25 A1+/? + 2n final states are so 
large that all calculated xn or xp (x * 1, 2) semi-inclusive 
cross sections significantly overestimate the present exper- 
imental measurements. In earlier work 15 on a semiempir- 
ical fragmentation code which used this same FSI model, 
it was noted that improved agreement between calcula- 
tions and all available experimental data were obtained if 
values of excitation energy were increased above those ob- 
tained from the methods of Ref. 11. In this work, we ob- 
served that treating £> si as a free parameter and increas- 
ing its value by 15 MeV reduced the cross sections for the 
25 Mg + 2 p + rt and 25 Al+/? + 2n channels — thereby 

improving the semi-inclusive cross-section predictions. 
Therefore, the final hadronic cross section for production 
of the type / isotope is given by 


m “n + z , 

(5) 

Zpf"Z— z , 

(6) 

A pf ^ A p m , 

(7) 


and ( m f ) denotes the usual binomial coefficient expression 
from probability theory. To complete the abrasion por- 
tion of the calculation, prefragment excitation energies 
must be estimated. We use 

£exc“£j+£*FSI » (8) 

where the surface energy term (£,) is calculated using the 
usual clean-cut abrasion formalism. 14 The frictional- 
spectator interaction (FSI) contribution (£>si) is estimat- 
ed using the methods of Ref. 11. To compute the proba- 
bility that p FSI’s have occurred for each abrasion of m 
nucleons, we use an extension of the Benesh, Cook, and 
Vary (BCV) prescription for estimating escape probabili- 
ties of abraded nucleons rather than the usual assumed 
value of one-half. 6, 11 Therefore, the abrasion cross section 
for a prefragment of isotopic species (Zpf,/1pf) which has 
undergone p FSI’s is given by 

Ot\>t(Zff,A Pr ,p)- (™](1 -P ac )f , (P a c) m - p 

*<Jibr(ZpF,Apf) , (9) 


&nuc(Zi ^ ^ Ct(j (.p ) Oabr (Zy , Aj ,/? ) , (10) 

J P ~ o 

where the summation over j accounts for the contributions 
to / from different prefragment species y, and the summa- 
tion over p accounts for the effects of the different excita- 
tion energies resulting from FSFs. 

Estimated exclusive cross sections obtained using the 
fragmentation model described herein are separately list- 
ed in Table I for each target. To compare our predictions 
with the semi-inclusive hadronic cross-section measure- 
ments (Fig. 4 of Ref. 1), we sum the exclusive channels 
listed in Table I for each of the relevant nucleon emission 
reactions. For example, the 1 p semi-inclusive calculation 
is the sum of the exclusive channel cross sections for the 
27 A1+/?, 26 A1 +/?+/t, 25 A1+/? + 2/i, and the 24 A1+/7 + 3 /i 
reactions. Similarly, the In semi-inclusive calculation is 
the sum of the 27 Si + n, 26 Al+p + n, 25 Mg + 2/? + n, and 
24 Na + 3/7 + n exclusive channels. The calculated results 
for xn (x-1, 2) and yp (y — 1, 2, 3) semi-inclusive cross 
sections are plotted in Fig. 1 along with the BNL experi- 
mental measurements from Fig. 4 of Ref. 1. Except for 
the 1/7 datum for the lead target, all calculated proton 
cross sections are in reasonably good agreement with the 
experimental data. Comparing the calculated and experi- 
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TABLE I. Exclusive channel hadronic dissociation cross- 
section calculations. 


Cross section (mb) with target nucleus 

Channel 

Aluminum 

Tin 

Lead 

27 Si+« 

99.1 

126.4 

134.4 

27 A1 +p 

99.1 

126.4 

134.4 

26 Si+2n 

17.9 

22.3 

23.9 

24 Al+p + rt- 

38.9 

48.2 

51.7 

26 Mg + 2p 

17.9 

22.3 

23.9 

25 A1+/? + 2/i 

1.5 

2.1 

2.3 

25 Mg + 2p + « 

13.8 

20.1 

21.8 

23 Na + 3p 

0.3 

0.4 

0.4 

24 AI +p + 3n 

0.1 

0.2 

0.2 

24 Mg + 2p + 2/t 

30.7 

44.2 

47.5 

24 Na + 3p + /t 

10.2 

14.5 

15.6 


mental neutron removal cross sections, however, we note 
that the agreement is not as good. There the calculations 
systematically overestimate the measurements by nearly 
50%. Since the experimental data were not corrected for 
detector acceptance limitations, 16 the observed trend for 
calculated cross sections to generally be larger than mea- 
sured ones is expected because the experimental data are 
likely to underestimate the actual cross sections by an as 
yet unknown amount. Resolution of these discrepancies is 
therefore hampered by the lack of exclusive channel mea- 
surements and detector acceptance corrections, which 
would enable the source(s) of any differences to be pin- 
pointed. 

For the 1 n removal calculations, the main contribution 
(nearly half) to these cross sections for each target arises 
from the 27 Si+/i exclusive channel when no FSI occur. 
Simple modifications to the current calculation model, 
such as using the Rasmussen 1 1 FSI escape probability 
(Pesc=0.5 ), would reduce the neutron cross-section 
differences; however, the calculated proton removal cross 
sections would also be reduced, destroying the agreement 
that presently exists between theory and experiment. A 
potential source for part of the difference between neutron 
and proton removal cross sections, not accounted for by 
the theory, is the difference in removal threshold energies. 
A proton, being less tightly bound, should have a larger 
removal cross section than a neutron. To test this hy- 
pothesis, we turn to the earlier fragmentation measure- 
ment of carbon and oxygen beams by Lindstrom et ai 17 
which provide a fairly complete data set. Correcting their 
measurements for EMD contributions using Ref. 18, we 
find that the exclusive 1 p removal channel ( ,5 N or ll B 
formation) is only 10%-20% larger than the exclusive 1 n 
removal channel ( l5 0 or !1 C). Adding the other In and 
1 p exclusive channels ( 14 N, l3 N, l3 C, ,0 C, etc.) to esti- 
mate experimental 1 n and 1 p inclusive cross sections 
yields much smaller differences between them — unlike the 
recent 28 Si measurements 1 where the 1 p semi-inclusive 
cross sections are substantially larger than the In cross 



Target mass number 

FIG. I. Hadronic dissociation cross sections vs target mass 
number. The experimental data point symbols include error 
bars; the theoretical calculation point symbols do not. I p is rep- 
resented by a solid square, In by an open square, 2 p by a solid 
circle, 2n by an open circle, and 3 p by a solid triangle. 


sections. From a binding energy point of view, this may 
result from the fact that the l |C- ll B binding-energy 
difference is smaller than that for 27 Si- 27 Al. A way to in- 
corporate proton-neutron binding-energy differences into 
the present model may be to use different nuclear distribu- 
tions for the proton and neutron densities. Such efforts 
are considered in Refs. 19 and 20. Recent work 21 22 has 
shown how the binding energy is directly influenced by the 
nuclear density. In principle, then, one could model the 
proton-neutron densities of 28 Si to fit the observed 
binding-energy differences. However, this particular 
method is beyond the scope of the present treatment. In- 
stead, possible changes to the calculated cross sections, re- 
sulting from neutron-proton density differences, were 
modeled by reducing the half-density radius of the 28 Si 
neutron distribution by 0.5 fm. The calculated neutron 
cross sections were reduced, as anticipated, but by only a 
few millibarns (less than 10 mb for all targets). These 
reductions were not large enough to account for the re- 
ported differences between measured semi-inclusive pro- 
ton and neutron removal channels. Clearly, exclusive 
channel experimental measurements for fc8 Si hadronic 
cross sections, which are presently being analyzed, 
would substantially aid efforts to resolve the differences 
between calculation and experimental measurements. 
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